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ABSTRACT 

We present a fast iterative FFT-based reconstruction algorithm that allows for non¬ 
parallel redshift-space distortions (RSD). We test our algorithm on both N-body dark 
matter simulations and mock distributions of galaxies designed to replicate galaxy 
survey conditions. We compare solenoidal and irrotational components of the redshift 
distortion and show that an approximation of this distortion leads to a better estimate 
of the real-space potential (and therefore faster convergence) than ignoring the RSD 
when estimating the displacement field. Our iterative reconstruction scheme converges 
in two iterations for the mock samples corresponding to BOSS CMASS DR11 when 
we start with an approximation of the RSD. The scheme takes six iterations when 
the initial estimate, measured from the redshift-space overdensity, has no RSD cor¬ 
rection. Slower convergence would be expected for surveys covering a larger angle on 
the sky. We show that this FFT based method provides a better estimate of the real 
space displacement field than a configuration space method that uses finite difference 
routines to compute the potential for the same grid resolution. Finally we show that 
a lognormal transform of the overdensity, used as a proxy for the linear overdensity, 
is beneficial in estimating the full displacement field from a dense sample of trac¬ 
ers. However the lognormal transform of the overdensity does not perform well when 
estimating the displacements from sparser simulations with a more realistic galaxy 
density. 
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1 INTRODUCTION 

Numerous upcoming galaxy surveys including WEAVE 
(Dalton et al. 2012), 4MOST (de Jong et al. 2012), DESI 


( Levi et al.|2013 ), eBOSS (Dawson et. al. in prep.), EUCLID 
( Laureijs et al.|20ll ) and WFIRST ( |Spergel et al.|2013 ) will 


measure the baryon acoustic feature to high precision at a 
range of redshifts. Their goal is to track the expansion his¬ 
tory of the Universe to sufficient precision to rigorously test 
cosmological models and decipher the mechanism behind the 
late-time accelerated expansion rate of the Universe. The 
Baryon Acoustic Oscillation (BAO) scale provides a stan¬ 
dard ruler in the distribution of mass and in turn galaxies, 
providing a tool to make such measurements. 

The BAO feature is created when sound waves prop¬ 
agating through the photon-baryon fluid stall as the pho¬ 
ton pressure is released during recombination (e.g. Meiksin, 


White & Peacock 1999). The stalled sound wave deposits 


baryonic perturbations in spherical overdense shells around 
dark matter perturbations. This generates a preferred scale 
in the matter overdensity distribution, the baryon acoustic 
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feature, large enough (« 120 h~ x Mpc) to remain detectable 
at low redshifts in the distribution of a large-volume sample 
of galaxies. Statistically the feature emerges in the 2-point 
correlation function as a bump or in the power spectrum as 
a series of peaks and troughs in the amplitude. 

Although the signal is robust, bulk flow motions and 
non-linear structure formation statistically blur the bound¬ 


aries of the shell increasingly to lower redshifts (Eisenstein, 
Seo & White 2007). The process of reconstruction recovers 


some of the lost precision due to these effects and is therefore 
crucial to fully exploit future galaxy survey data. 

The simple reconstruction method of Eisenstein et al.| 


Anderson et al. 2014a 

Tojeiro et al. 2014||Ross et al. 

2014 

Anderson et al. 2014 

o\ Kazin et al. 2014| Padmanabhan| 


the baryon acoustic signal. The method relies on the abil¬ 
ity to recover the displacements of overdensities from their 
initial positions and reverse this motion. Practically this is 
achieved by moving galaxies and a set of random particles 
that trace the geometry of the survey back along these vec¬ 
tors. Statistically this sharpens the measured signal and has 
been shown to remove a slight bias in the distance measure- 
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ment induced by non-linear structure formation (Eisenstein, 
Seo fe White|2007| [Mehta et al.|2011 ). 


It is straightforward to compute an estimate of the dis¬ 
placement field from the smoothed real space overdensity 


using the Zel’dovich approximation (Zel’dovich 1970) and 


Fast Fourier Transforms (FFTs). However, complications 
arise due to RSD which distort the measured overdensity 
from the true overdensity. One approach to include RSD is 
to use finite difference methods to linearly model both real 
space potential and RSD (Padmanabh an et al.| [2012). The 
real space potential is recovered at each location from the 
overdensity field measured in redshift-space using a linear 
equation solver. The displacement field is computed from 
the potential using finite differences. The accuracy of this 
method is limited by the grid resolution adopted. Further¬ 
more, solving for the potential on the grid using a linear 
equation solver is computationally expensive. 


Although one cannot directly compute the real space 
displacement field from the redshift-space overdensity using 
FFTs because of the varying line of sight, |Burden et al.| 
(2014) made an approximation that allows partial correction 


of the redshift-space component. The resulting displacement 
vectors were shown to be well matched to those computed 
in configuration space. We now show that we can improve 
upon this by developing an iterative scheme to correct for 
RSD. 


Both FFT and finite difference methods require us to 
use the overdensity field to estimate the displacements. 


Neyrinck, Szapudi & Szalay (2009) show that the lognor¬ 


mal transform of the overdensity field measured in N-body 
simulations reduces the non-linear component of the field on 
quasi-linear scales thereby recovering more linear informa¬ 
tion in the power spectrum. The lognormal transform of the 
field is more Gaussian, thus off diagonal elements of the co- 
variance matrix are reduced. In Seo et al. (2011) the authors 
apply the lognormal transform to the weak-lensing conver¬ 
gence field and find similar results. The nonlinear informa¬ 
tion in the convergence field is reduced and the information 
content in the covariance matrix enhanced. Furthermore it 
has been shown (Falck et al.|[2Q12| ) that the divergence of 
the true displacement field (using N-body simulations) is 
better replicated under certain conditions by a logarithmic 
transform of the overdensity. We now test whether this im¬ 
provement also helps with reconstruction. 


In Section [2] we outline two basic methods used to com¬ 
pute the displacement vectors in the reconstruction process. 
In Section [3] we describe the data sets used to test our algo¬ 
rithms. In Section^] we test how well the the displacements 
field computed from the methods in Section [2] match the real 
space displacement field. In Section [5] we present our itera¬ 
tive reconstruction algorithm. In Section [6] we test the effect 
of using a lognormal transform of the overdensity to compute 
the displacement vectors in real space. Finally in Section [7] 
we conclude by combining the results of our tests to offer 
our most effective method of recovering the real space first 
order Lagrangian displacement field from an evolved galaxy 
over density field measured in redshift-space. 


2 RECONSTRUCTION METHODS 


The reconstruction process requires computing the La¬ 
grangian displacement field from the evolved galaxy over¬ 
density distribution measured in redshift-space. In this sec¬ 
tion we outline two methods to do this. 

A Lagrangian co-ordinate system tracks the motion of a 
fluid element through space and time where the Lagrangian 
displacement vector \l/ (q, t) maps a fluid at some initial po¬ 
sition q to its Eulerian position some time later x(q, t), 

x(q,t) = q+^(q,i). (1) 


The overdensity field is defined as <5(x) = (p(x) — p)/p and 
relates the density, p, at location x to the expected average 
density, p. 

The Zel’dovich approximation (Z el’dovich| [l970) is the 
first order term in a perturbative expansion of \I/(q, t) in 
Lagrangian perturbation theory (Zel’dovich||1970 |Buchert 
19931 P~994l |Bouchet et aL]|1995a ~ |Catelan||1995| Taylor k, 


Hamilton|1996 ), 


*(q) 


(i) _ 


/ 
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The successive perturbative terms T' fri) are functions of 
powers of the linear overdenisty perturbation. We use the 
Zel’dovich approximation throughout this work, and there¬ 
fore drop the superscript. 

To extract the displacement vector field from a redshift- 
space galaxy overdensity field, one must solve the equation 
(Nusser & D avis| 1994) 

_ c 

V-V + fV-(V-r)r=—. (3) 


This can be solved on a grid in configuration space, comput¬ 
ing gradients using finite difference approximations. How¬ 
ever, solving for a value of \l/ at each grid point is com¬ 
putationally expensive and the accuracy of the method is 
dependent on the resolution of the grid. Alternatively, [Bur] 
den et al. (2014) made the assumption that (\P • ?)? is an 


irrotational field and use FFTs to compute the value of \l/ 
increasing the speed of the computation. We now outline 
these two methods. 


2.1 Finite difference approximations 


The configuration space method solves for the potential, 0. 
Assuming \l/ is an irrotational field, as expected on linear 
scales, this is related to the first order Lagrangian displace¬ 
ment as \l/ = —V0 and Eq. [3] can be written as 

V 2 ^ + /V-(V^)f = (4) 


We follow Padmanabhan et al. ([2012) and solve this on a 
grid using finite differences to approximate the derivatives. 
The potential at each grid point is expressed as a function of 
the potential at the surrounding grid points. The Laplacian 
of the potential at a grid point can be approximated as a 
function of the potential at the 6 nearest grid points and the 
central point 


V 0ooo 


1 

g 2 


<fiijk — 60000 , 
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(5) 
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where the sum over A is the sum over the 6 adjacent grid 
points and g is the spacing between grid points. The second 
part of Eq. [4] can be written as 

/v • (V0 r ) f = / (? • V (V0 r ) + V0 r (V • f)), (6) 


which can be approximated as 
.0000 


r 


- 2 ^ + £/^ + £(-1)'/ 


0^' 




2p 2 


(7) 

where B is the set of points ijk such that 2 of the indices 
are zero and the other is ±1, is the cartesian position of 
the non-zero index and r is the distance to the central grid 
point. C is the set of points where two of the indices are ±1 
and the third is zero. When the two indices are the same, 
p = 0, when they are different p = 1. x\ and Xj are the 
cartesian positions of the non-zero indices. 

This can be arranged as a linear system of equations 
such that A0 = S g /b, where A is a matrix describing the 
dependence of the potential on its surroundings. The S that 
we input here is the smoothed overdensity field measured in 
redshift-space. We solve for the potential using the GMRES 


in the PETSc package (Balay et al. 2014 2013). Once we 


have a model for the potential, finite differences are used 
again to calculate the displacements at each grid point from 
the potential. 


improved by accounting for the solenoidal component. In 
Section 5 we show that the approximations do not change 
the outcome of our new iterative process designed to solve 
Eq. [3] numerically, although they do speed up convergence. 


2.3 Displacement field nomenclature 

As we use a number of estimates of the displacement field 
computed using a variety of methods, we introduce our no¬ 
tation here; 

• ^true, the true displacement field measured from final 
minus initial particle positions (N-body only). 

• ^red, fd displacements computed from the redshift- 
space overdensity in configuration space using the finite dif¬ 
ference (FD) method. 

• ^real, fd displacements computed from the real space 
overdensity in configuration space. 

• real, fft displacements computed using FFTs and the 
real space overdensity. 

• red, fft displacements computed using FFTs and the 
redshift-space overdensity (with no RSD correction). 

• ^ffta displacements computed using FFTs and the 
redshift-space overdensity where an approximate correction 
assuming that (\k • r)r is irrotational has been applied. 


2.2 Fast Fourier Transform method with 
redshift-space approximation 

While is expected to be irrotational, (\k • ?)? is not. This 
means we cannot easily solve Eq. [3] using Fourier methods. 
To see this we decompose (\k • ?)?, using Helmholtz’s Theo¬ 
rem, into the gradient of a scalar potential field and the curl 
of a vector field 


(^•f)f = VA + VxB, (8) 


where A is a scalar potential held and and B is a vector 
potential held. We refer to VA as the irrotational component 
and V x B as the soleniodal component. Substituting Eq. [8] 
into Eq. [3] and writing as the gradient of a potential as in 
Eq. [4] we get 

V(0 + /A) = -VV- 2 ^. (9) 


This is exact and without approximation, but does not allow 
us to calculate directly using FFTs. A further assumption. 
Burden et al. (2014) made the approximation that (\k • ?)? 


is irrotational, that is 


($-f)f«VA, (10) 

making Eq. [9] 

'*' + /('*'• f)f=-VV _2 |. (11) 

The RHS can be computed in Fourier space using FFTs, 
and the solution to this equation is 


^ = IFFT 


(k) 

k 2 b 


f 

1 + / 


jlFFT 


ik£ (k) 
k 2 b 



where IFFT is the inverse Fourier Transform. 

In Section we show that this approximation can be 


3 CATALOGUES 

In this section we describe the two sets of simulations used 
to test our algorithms. 


3.1 N-body simulations 

We use a TreePM N-body simulation box of 1024 3 dark 
matter particles generated with GADGET-2 (Springel;2005 


Springel, Yoshida & White 2001) with initial positions at 


z = 49 and hnal positions and velocities at z = 0. The box 
has side length 768 h~ x Mpc comoving coordinates. The fidu¬ 
cial cosmology of the simulation has = 0.317, h = 0.670, 
Q b h 2 = 0.022. 

The full sample is randomly subsampled to 20(128 3 ), 
5(128 3 ), 128 3 and 0.25(128 3 ) particles to create simulations 
with varying density. All of these samples are derived from 
a single realisation. For brevity these are named the 20t, 5t, 
It and 0p25t samples respectively. In order to investigate 
redshift-space effects we locate an observer at z = 0, and 
modify the radial positions by (v • f )/Ho where v is the final 
velocity in units kms -1 and the Hubble constant is Ho is 100 
kms -1 h Mpc -1 and r is the unit vector pointing away from 
the origin of the box. This gives a sample with strong wide- 
angle effects, a very wide sky area (1/8 of the sky) near to 
the observer. 

The particle overdensity is measured as described in 
Section [ 2 ] where p(x) is the density of the particle distribu¬ 
tion at grid point x and p is the average particle density 
per grid point. We use a NGP binning scheme on a 128 3 
mesh. We do not investigate alternative binning methods as 
the grid points are smoothed post binning mitigating any 
effects of different binning schemes. 

Gaussian smoothing filters are applied in Fourier space 
as a multiplicative term S(k) = exp— ik 2 R 2 /2 where R is 
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Figure 1. The plots show a comparison of irrotational and solenoidal components of (SI* • x)x projected in cartesian directions. The 
redshift space distortions are modelled as plane-parallel and projected along the x direction. The components projected along x (left) 
are positively correlated and the amplitudes are higher than in the other two directions. The average ratio of these vectors is 0.75, this is 
shown by the black line. The sum of the two components projected in the y (centre) and z (left) is expected to be zero, the components 
are negatively correlated in these two directions. The amplitudes of these components is small. 


the smoothing length. Smoothing removes non-linearities in 
the overdensity, reduces shot noise and sets the scale of the 
overdensities to be moved during the reconstruction process. 
The smoothed overdensity is used as a proxy for the linear 
overdensity to calculate the first order Lagrangian displace¬ 
ment field in Eq. [2] 


3.2 Mock survey catalogues 


To conduct our tests on galaxies with realistic masks, we 


make use of the publicly available PTHalo (Manera et al. 
[20l3| > mock samples created to match Data Release 11 
(DR11) of the Baryon Oscillation Spectroscopic Survey 
(BOSS) galaxy samples ( Dawson et al.||2013 ). 

We use the CMASS Northern sample spanning a red- 
shift range of 0.43 ^ z ^ 0.70 and covering an effective 
area of 6,308 square degrees. The PTHalo method initially 
creates a matter field based on second order Lagrangian Per¬ 
turbation Theory (LPT). Halos are located with a Friends 
of Friends (FoF) algorithm, and halo masses calibrated to 
N-body simulations. The halos are populated with galax¬ 
ies using a Halo Occupation Distribution (HOD) calibrated 
by the observed galaxy samples on small scales between 
30/i -1 Mpc and 80/i _1 Mpc using the 2-point statistics of 
the true galaxy sample. Redshift-space distortions are added 
to the mock galaxy distribution by modifying their redshifts 
according to the second order LPT peculiar velocity field in 
the radial direction. The matter field is created in a single 
time-slice, rather than in a light cone, thus the growth rate 
and RSD signal are constant throughout the sample. The 
effective redshift of the survey is z = 0.57. 

The mock catalogues are sampled to match the angular 
mask and redshift cuts of the survey data. To replicate ob¬ 
servational complications, galaxies are subsampled to mimic 
missing galaxies caused by redshift failure, and close pairs as 


simultaneous spectroscopic observations are limited to ob¬ 
jects separated by greater than 62". 

The fiducial cosmology of the mock sample has — 
0.274, h = 0.70, Q b h 2 = 0.022, n s = 0.95 and cr 8 = 0.80. We 
assume a local deterministic galaxy bias such that 5 g = bS, 
where b — 1.87 is the galaxy bias estimated empirically from 
the real survey data. 

As the distribution of mock galaxies includes sampling 
and mask effects of survey data we compute the overden¬ 
sity using a random catalogue of particles Poisson sampled 
within the survey mask. To minimise shot noise the ran¬ 
dom catalogue contains 100 times more data points than 
the mock catalogue. 

Although we do not compute the power spectrum in 
this work, we weight mock galaxies using the FKP weight¬ 
ing scheme in Feldman, Kaiser & Peacock (1994), designed 
to optimally recover the overdensity field given a varying 
density sample. We therefore apply a weight to each galaxy 


W — W F KP (w cp + W red - 1 ) , 


(13) 


where w cp and w re d correct for the close-pairs and redshift 
failures respectively (see Ande rson et al.||2014b for further 
details), and wfkp is the FKP weight 


wfkp = 


1 

1 + n (z) Po ’ 


(14) 


with fixed expected power spectrum Po = 20, 000h _3 Mpc 3 , 
and average galaxy density n(z). 

The galaxies are set in a box of length 3500/i _1 Mpc 
and the overdensity is computed using the NGP scheme on 
a 512 3 mesh. The mock catalogue is sparse, therefore we 
smooth the galaxy and random catalogue density prior to 
computing the overdensity. To prevent spurious spikes in 
our overdensity post smoothing, we multiply the smoothed 


© RAS, MNRAS 000, [lp?? 


















Reconstruction in Fourier space 5 


overdensity by a binary mask that sets all grid points that 
fall outside of the survey region to zero. The first order dis¬ 
placement held is computed as before. 


4 THE EFFECTS OF AN IRROTATIONAL 
APPROXIMATION 


In this section we investigate the RSD term in Eq. @ 0-r)r, 
and in particular consider the approximation that it is irro- 
tational. We show that this approximation can be improved 
beyond that of Burden et al. (2014). 


4.1 The plane-parallel framework 

Although 4/ is irrotational ( Bouchet et al.| 1995b), (4> • r)r 
is not and the true vector held has both irrotational and 
solenoidal components as shown in Eq. [8] To estimate the 
amplitude of these components we assume a distance ob¬ 
server plane parallel approximation of the RSD where r —> 
x. Taking the curl of (4/ • x)x gives 


V x (4> • x)x = V x (V x B). 


(15) 


Therefore in Fourier space the solenoidal components of the 
redshift-space distortion of the displacement held are 

VxB = ± (k 2 v + kl) ^k x -^^k y -^^k z , (16) 
and the irrotational components are 

+ + (17) 

where we have made the abbreviation = 4> • x. We see 
that V x B is non-zero and has components of the same 
order of magnitude as VA. As expected, off-axis components 
cancel in the sum. 

To empirically measure the amplitude of these compo¬ 
nents we compute the three displacement held matrices that 
comprise 4> tr ue at each grid point in the N-body catalogue. 
We smooth these using a Gaussian hlter with smoothing 
scale R — 10 h _1 Mpc[jWe then compute the soleniodal and 
irrotational components of (4/ • x)x in Fourier space using 
FFTs as in Eq’s. [16] and [lT| The measurements are con¬ 
verted back to configuration space for comparison. Fig. [I] 
shows scatter plots comparing the amplitudes of the vector 
helds in cartesian directions. A histogram of the ratios of 
irrotational to solenoid components is shown in Fig. [2] 

As expected from Eq’s. [l6| and [l7| (VA)-y and (VxB)- 
y are negatively correlated and (VA) z and (VxB)-z are 
negatively correlated as the RSD only lie along the x-axis. 
The amplitude of these terms is small compared with that 
along the x-axis. 

The Helmholtz components projected along x are pos¬ 
itively correlated and have a larger amplitude. In the plane 
parallel approximation the irrotational and solenoidal com¬ 
ponents projected along x have a mean ratio of 0.75. From 
this empirical result we make a new approximation that 


1 This smoothing length chosen is shown to perform best when 
reconstructing the linear BAO signal extracted from spherically 


averaged power spectrum measurements (Burden et al.|2014). 



Figure 2. A histogram of the ratios of irrotational to solenoid 
components of (4/ • x)x measured in our N-body simulation. The 
black line is the ratio projected along the x axis with mean value 
of 0.75 and standard deviation, cr = 2.50. The blue line shows the 
values projected along y with mean of -0.75 and cr = 2.25, and 
the red line shows the same projected along z with mean of -0.85 
and cr =1.95. The differences between the blue and red curve are 
due to noise as only one simulation is used in this test. 


VA • x « 0.75(V x B) • x. With this approximation, as¬ 
suming that the geometry of the survey allows the plane 
parallel approximation, we have that VA « (3/7) (4/ • r)r so 
that Eq. [TT] becomes 


f+^(f-f)f = -VV 2 ^, (18) 

rather than Eq. ED This has the solution 


^ = IFFT 


ik5g (k) ‘ 

f IFFT 

*k5 g (k)" 

k 2 b 

7 + 3/\ IFFT 

k 2 b 


r. 


(19) 


This suggests that our initial approximation in Eq. [12] over 
compensates for the RSD component of the displacement 
held. This approximation is derived empirically from one 
N-body realisation. To use this method for data analysis, 
tests are required over a large number of realisations with 
varying cosmological parameters. In Section [5] we show that 
the approximations are not important for the end result if 
we adopt an iterative approach. They only affect speed of 
convergence. 


4.2 More general geometries 

We compare the real displacement held projected in the ra¬ 
dial direction, (4/reai, fft •?) to the displacements computed 
using the non-corrected held (4/ re d, fft • f). Then we com¬ 
pare the real displacement held projected in the radial di¬ 
rection to the displacements computed using the approxi¬ 
mations in Eqs. |12| and |19| projected in the radial direction, 
(4>ffta • r), using the N-body catalogue. The comparisons 
are quantihed by taking the mean ratio of these values at 
each grid point. The mean ratios and standard deviations 
are shown in the top part of Table [l] To prevent true dis¬ 
placements close to zero skewing the distribution, we cut out 
outlying ratios that are above 3 or below -1. The percent¬ 
age of outliers in each computation is shown in the far right 
column of the table. 

With no correction, the mean ratio of the amplitude of 
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/ ^real, FFT ‘ r \ 


outliers [%] 

\ / 


N-body 

^red, FFT _ 

0.79 

0.28 

2 

T r FFTA(Eq. 12) 

1.13 

0.34 

3 

^FFTA(Eq. 19) 

0.96 

0.31 

2 

Galaxy mock 

'I'red, FFT _ 

0.80 

0.42 

4 

\I>FFTA(Eq. [19]) 

1.04 

0.47 

5 

^red,FD 

1.12 

0.48 

6 

^real,FD 

0.97 

0.52 

7 


\I> 

/ ^real, FFT ' U \ 
\ * * T ± / 

c t outliers [%] 

\I>FFTA(Eq. 19 

) 1.05 

0.53 2 


Table 1 . Mean ratio and standard deviation of the the real space 
displacement field to the displacement field computed using dif¬ 
ferent estimations (see subsection |2.3| for notation) projected in 
the radial direction for both N-body simulations (top) and mock 
galaxies (centre). Outliers, defined as ratio values above 3 or less 
than -1 have been cut from the analysis. The average values of 
the mean ratio and standard deviation between the estimated dis¬ 
placement field and real space displacement field projected in the 
off-radial directions (rj_) is shown in the bottom row. 


the radial component of the real displacement field to that 
of the redshift-space displacement field is 0.79 with a stan¬ 
dard deviation of 0.28. The enhanced clustering in the radial 
direction in redshift-space leads to an overestimation of the 
amplitude compared to the real space displacement vectors. 
Using the approximation in Eq. ^]the mean ratio between 
the amplitude of the real space field projected in the radial 
direction and the corrected displacement field is 1.13 with 
a standard deviation of 0.34. As predicted the displacement 
field in the radial direction is overcorrected. Using the ap¬ 
proximation in Eq. [l9]the mean ratio between the amplitude 
of the real space field projected in the radial direction and 
the corrected displacement field is 0.96 with a standard de¬ 
viation of 0.31. The amplitudes of the displacements are on 
average 5% larger than those of the real space displacement. 
The approximation in Eq. [19] therefore brings the estimate 
of the displacement field closest to the real space displace¬ 
ment field even though the N-body simulation does not have 
a plane-parallel RSD geometry. 

Next we use the mock galaxies to compare the real dis¬ 
placement field projected in the radial direction, (\Eheai, fft * 
f) to the non-corrected displacement field, (SEhed, fft • r), 
and the displacement field computed using the approxima¬ 
tion in Eq.[l9] (^ffta • f) projected in the radial direction. 
Following this we compute the displacement field using the 
finite difference method outlined in subsection |2.1| from the 
redshift-space overdensity. In order to compare the results 
of both methods in real and redshift-space we also compute 
the real space displacement field using finite differences. 

With no correction, the mean ratio of the amplitude 
of the radial component of the real displacement field to 
that of the redshift-space displacement field is 0.80 with a 
standard deviation of 0.42. The ratio shows that (as seen in 
the N-body simulation), the displacements computed from 
the redshift-space data are larger, on average, than the 


real space displacement vectors. Using the approximation 
in Eq. ^]the mean ratio between the amplitude of the real 
space field projected in the radial direction and the corrected 
displacement field is 1.04 with a standard deviation of 0.47. 
The approximation brings the estimate of the displacement 
field closer to the real space displacement field. The ampli¬ 
tude of the radial displacements are smaller by 4% than the 
real space values. 

The mean ratio of the real space displacement fields 
computed using FFTs and finite differences projected in the 
radial directiorois 0.97 with a standard deviation of 0.52. 
Therefore there is only a small difference between methods 
applied to real space data as expected. The mean ratio of 
the amplitude of the real space displacement field and the 
redshift-space displacement field both computed in configu¬ 
ration pace and projected in the radial direction is 1.12 with 
a standard deviation of 0.48. Therefore the finite difference 
method overcorrects the redshift-space displacements in the 
radial direction. The accuracy of finite difference scheme is 
dependent on the resolution of the grid used. Finer grid res¬ 
olution would improve the accuracy of the computations, 
however, the procedure would be more computationally ex¬ 
pensive. The average mean ratios and standard deviations 
are shown for all of these results in the second part of Ta¬ 
ble □ 


Finally, to check that the approximation in Eq. |19| does 
not induce distortions in the off-radial directions (?_l), we 
compute the ratios (SEheai, fft • £_l) to (^ffta • ?_i_). The 
mean ratio is 1.05 with a standard deviation of 0.53. The 
values are the same for both off-radial directions. We con¬ 
clude that any distortions cast by the approximation in the 
non-radial directions are negligible. The values are tabulated 
at the bottom of Table |T] 

In this section we have shown that the approximation 
(^ r)r is irrotational as adopted by |Burden et al.| p014) 
overcorrects for the RSD component of the displacement 
field. However it is closer to the real space displacement 
field than the non-corrected field. Assuming a plane parallel 
approximation of RSDs we empirically derived the relation¬ 
ship between the irrotational and solenoidal components and 
used this to improve our estimate the real space displace¬ 
ment field. Our updated approximation is shown in Eq. |19| 
We show that this approximation, computed using FFTs, 
is closer to the real space displacement field than one com¬ 
puted in configuration space which overcorrects the RSD 
component of the displacement field. 


5 AN ITERATIVE RECONSTRUCTION 
SCHEME 

In this section we propose an iterative method of estimating 
the real space overdensity from the redshift-space overden¬ 
sity field. Note that this iterative method is a precursor to 
the actual reconstruction process. We apply the iterative 
method to both the N-body and mock catalogues. The ra¬ 
tio of the amplitude of the real space overdensity field to 
the amplitude of the displacement field computed at each 


2 The field is isotropic but we chose the radial direction for con¬ 
sistency. 
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Figure 3. Iterative convergence of the ratio of (SE'real, fft • r) to (^est • r) measured from the N-body simulation (left) and the mock 
galaxy catalogue (right). Starting with the non-corrected field (red dashed line), the N-body ratios at each iteration oscillate and converge 
after 4 iterations. Starting with the corrected field (blue full line), the values converge after 2 iterations. Starting with the non-corrected 
field (red dashed line), the mock catalogue ratio at each iteration oscillates about 1 and converges after 6 iterations. Starting with the 
corrected field (blue full line), the mock catalogue values converge at 2 iterations. The criterion to define convergence is stated in the 
main text.The same ratio from the finite difference computation is 1.12 (grey point). 


iteration in the radial direction is compared. This iterative 
method allows us to avoid the approximations discussed in 
the previous section. 


5.1 Method 

As seen in Section [4] we can estimate \l/ using Eq. |19[ or we 
could have ignored the RSD completely and our estimate 
would be \l/ = IFFT [ik$ gj red /k 2 bj . In terms of the poten¬ 
tial, our estimate is 0 es t, 1 = IFFT £ gj re d/k 2 6]. From Eq.| 
we expect our estimate to be related to the potential 0 by 

0 est, 1=0 + ( 20 ) 

To simplify the algebra, we have defined an operator 

G{4>) = V _2 V • (V0-f)f. ( 21 ) 

The estimate of the potential can be used to remove the 
redshift component of the displacement field so that our first 
estimate of the real space overdensity is 

<5g ’ re b a1 ’ 1 = + /V • (V (test, 1 • f)f, ( 22 ) 

which, substituting for 0 es t, i becomes 

= _v 2 0+/ 2 V • (VG(4>) -f)r. (23) 

From this we can make a second estimate of the potential 

^est .2 = -V- 2 +^l=0-/ 2 G(G(0)). (24) 

As before, we use it and Eq. @ to remove the redshift compo¬ 
nent of the displacement field giving us our second estimate 
of the real space overdensity 

5jtre b a1 ’ 2 = +P 1 + /V • (V <pest, 2 • f)f. (25) 

Or equivalently 

<5g ’ raa1 ’ 2 = -v 2 0 - / 3 V • (VG(G( 0 )) • f)f, ( 26 ) 

and we see our recursive iterative method emerge. In general 
for each iteration n, the corrected overdensity is related to 


the true galaxy over density by 

<W, n + (-/)"+i V ■ (VG (n) (0) • r)r, (27) 

where the subscript n in the Gr( n )(0) term describes the 
number of times the function G(0) is applied to itself. 
We then go on to use £9,real, n to estimate \I/fft, n — 

W ($g, real, n/b'). 

For a plane parallel approximation in Fourier space 
where f — y x the n th iteration of the displacement field can 
be written as 






FFT, n — 


'g, real 


bk 2 


1 + (-/) (n+1) (+) 


2(n+l) 


(28) 


For the iterative process to converge, the second term on the 
RHS of Eq. [28] must decrease in magnitude as n increases. 
As k x ^ k it will be reduced in amplitude at each iteration 
provided / < 1. Similar convergence is expected for more 
general geometries given that the relative amplitudes of the 
terms do not change. 


5.2 Results 

We consider two starting points for the iterative procedure. 
The first is the displacement field computed directly from 
the redshift-space coordinate using FFTs, \l/red, fft- The 
second is the same field modified by the approximate cor¬ 
rection given in Eq. [19] We call this \I/ffta. The iterated 
versions of these fields are denoted \I7fft, n and \I/ffta, n 
respectively where n is the iteration number. From Eq. [27] 
we see the ratio of the displacement field measured from the 
real space overdensity projected in the radial direction to 
that of the displacement field computed from the redshift- 
space overdensity projected along the radial direction will 
be 

real, FFT • f _ _ ^real, FFT ‘ 

^FFT,n • f ^real, FFT * r + ( —/)+ + 1 )G^ • ? ’ 

where, if the term 

G* = -VV 2 [V • (VG (n) (0) • f)f] (30) 
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^est 

n 

/ ^real, FFT 

\ 4>est • r 

-> • 

outliers [%] 

^FFT, n 

0 

0.79 

0.28 

2 


1 

1.07 

0.50 

7 


2 

0.97 

0.32 

3 


3 

1.01 

0.40 

4 


4 

0.99 

0.36 

3 


5 

1.00 

0.37 

4 

\I>FFTA, n 

0 

0.96 

0.31 

2 


1 

1.02 

0.44 

5 


2 

0.99 

0.35 

3 


3 

1.00 

0.39 

4 


4 

1.00 

0.37 

4 


5 

1.00 

0.37 

4 


Table 2. Iterative results for the N-body simulation. There are 
20(128 3 ) particles in the sample which has been placed on a 128 3 
grid and smoothed with R=10 h -1 Mpc. The mean ratio and stan¬ 
dard deviation of the ratio of the real space displacement field 
projected in the radial direction to that computed at each iter¬ 
ation are shown. The two starting points are the displacement 
field measured directly in redshift space, ^fft, (o) an d the ap¬ 
proximation in Eq. fl 9 l 'Fffta, (o)- We define convergence as the 
iteration to which subsequent values of the mean differ by less 
than 1%. Starting with no approximation the iterative procedure 
converges in 4 steps. Starting with the approximation the pro¬ 
cedure converges in 2 steps. Outliers are defined as ratio values 
greater than 3 or less than -1. 


^est 

n 

/ ^real, FFT 
\ ^est • r 

-> * 

outliers 

^FFT, n 

0 

0.80 

0.42 

4 


1 

1.10 

0.51 

7 


2 

0.96 

0.43 

4 


3 

1.03 

0.47 

6 


4 

0.98 

0.43 

5 


5 

1.02 

0.46 

5 


6 

0.99 

0.44 

5 

^FFTA, n 

0 

1.04 

0.47 

5 


1 

1.02 

0.46 

5 


2 

1.00 

0.43 

5 


3 

1.01 

0.45 

5 


4 

1.00 

0.44 

5 


5 

1.01 

0.45 

5 


6 

1.00 

0.44 

5 


Table 3. Iterative results for the mock catalogue. There are 
601921 galaxies in the sample which has been placed on a 512 3 
grid and smoothed with R=10h -1 Mpc. Shown is the mean ratio 
and standard deviation of the ratio of the real space displace¬ 
ment field projected in the radial direction to that computed at 
each iteration projected in the radial direction. The two starting 
points are the displacement field measured directly in redshift 
space, ^fft, (o) an d the approximation in Eq. 1191 ^ffta, (o)- 
We define convergence as the iteration to which subsequent values 
of the mean differ by less than 1%. Starting with no approxima¬ 
tion the iterative procedure converges in 6 steps. Starting with 
the approximation the procedure converges in 2 steps. Outliers 
are defined as ratio values greater than 3 or less than -1. 


decreases in magnitude with each iteration, the ratio will os¬ 
cillate around one and converge provided we use the correct 
growth value in the procedure. 

In Fig. [3] we show how the mean of the ratio in Eq. [29] 
computed at each particle/galaxy position changes for suc¬ 
cessive iterations in both N-body and mock catalogues. For 
comparison we also show the ratio computed from the finite 
difference algorithm applied to the mock catalogue. The val¬ 
ues are provided in Table [2] (N-body) and Table [3] (mock). 

As predicted, the mean values of the ratios (\tReai, fft • 
r) to (\l/FFT,n • r) oscillate around one and converge. Note 
that we compare smoothed fields and the smoothing is a 
convolution that does not cancel out in the ratio. We define 


6 LOGNORMAL TRANSFORM OF S 


It has been shown (Falck ~ al.|2012 ) that under certain con¬ 
ditions, a logarithmic transform of the overdensity provides 
a better estimate of the divergence of the displacement field 
in N-body simulations. In ( Falck et al.| 2012), the smooth¬ 
ing of the overdensity field is achieved using a CIC binning 
assignment, however we use a Gaussian smoothing kernel in 
this work, the goal of which is to explore the effect of sample 
density of the lognormal results. It should be noted that our 
densest sample is approximately 20 times less dense than the 
simulations used in their work. In this section we consider 


convergence as the point at which further iterations do not 
change the mean value of the ratio by more than 1%. There 
is no further improvement after iteration number 4 for the 
N-body simulation and iteration number 6 for the mock data 
sets. In the N-body catalogues the mean ratio converges to 
0.99 with a standard deviation of 0.36. In the mock galaxies 
the mean ratio converges to 0.99 with a standard deviation 
of 0.44. Starting the iteration with the corrected field, 2 it¬ 
erations are needed for the ratio to converge to 0.99 with 
a standard deviation of 0.35 in the N-body simulation. Us¬ 
ing the galaxy mock sample 2 iterations are needed for the 
mean ratio to converge to 1.00 with a standard deviation of 
0.43. As in Table [l] the column on the far right shows the 
percentage of outliers that have a ratio greater than 3 or less 
than -1 and have not been included in the computations. 

Fig. [4] compares the amplitude of \EReai, fft • ? to 
^ red, FFT * ?, ffta • f and \I/ffta,2 • f computed from 
the mock catalogues at each galaxy position. It is seen that 
the tightest correction with the real space displacement field 
comes from the \I/ffta,2 • r field. 


if such a transform can help with our reconstruction algo¬ 
rithm using different density N-body samples. We first test 
the transform in the calculation of V • T r . We then test if the 
lognormal transform combined with the Zel’dovich approxi¬ 
mation improves our estimate of \l/. All of the comparisons 
are performed in real space. 

The Newtonian continuity equation relates the rate of 
change of a density perturbation to its peculiar velocity v = 
ax. It can be expressed (retaining the complete (1 + 5) term 
in the time derivative ( |Pee b ies||1980| )as 


0(1 + 6 ) 
dt 


+ -V • (1 + S)v = 0. 


(31) 


This equation can be linearised about the quantity S or (1 + 
5). The respective linear continuity equations read 


OS 1 

— + - V-v = 0, 
dt a 


(32) 


1 0(1 + 6 ) 

(1 + S) dt 


+ iv-v = 0, 

a 


(33) 
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no. gals 


no. gals 


no. gals 



Figure 4. From left to right the plots show a comparison of the real space displacement vectors projected in the radial direction to; 
the field computed directly from the redshift-space overdensity ^red, FFTl the field computed with the approximation given in Eq. |19| 
'Fffta; and the field computed by iterating twice starting with the same approximation, \I>ffta,2- The redshift-space field overestimates 
the displacement field, the approximation moves the two distributions closer together and the iterated field shows the tightest correlation 
with the real space field. The 1:1 values are described by the black lines. 


Sample 

Mpc] 

°A 

&B 

20t 

10 

0.57 

0.43 

5t 

10 

0.61 

0.49 

It 

10 

0.79 

0.70 

0p25t 

10 

1.27 

1.21 

0p065t 

10 

2.26 

2.48 

20t 

5 

1.01 

0.74 

5t 

5 

1.05 

0.82 

It 

5 

1.28 

1.16 

0p25t 

5 

1.89 

2.15 

0p065t 

5 

3.28 

5.66 


Table 4. The standard deviation computed form the distribu¬ 
tions 4>true — 4>(5) denoted a a and 4>true — 4>(ln[l + 8] + C) 
denoted as- As the density of the sample is decreased, the im¬ 
provement seen by using the the lognormal transform of the over¬ 
density to compute the displacement field decreases. For both the 
5 h _1 Mpc and 10 h~ x Mpc smoothing scale at the lowest density 
sampled, the lognormal transform performs worse than the over¬ 
density when used in the ZePdovich approximation. 


and have solutions 


< 

ii 

i 

(34) 

(1) = — ln(l + <5 (1) ) + C. 

(35) 


V ^ 


For each N-body density subsample, the overdensity is 
computed and smoothed with R — 5/i -1 Mpc and R — 
10/i -1 Mpc. The lognormal transform of the overdensity 
field is computed, ln(l + S a (5) + C using the smoothed 
8 values. Following Falck et al. (2012) we choose C = 


— (ln(l + S * 6)). The displacements, \F t rue computed us¬ 
ing the 20t sample are smoothed with the same values of 
R. The divergence of the displacement field is computed in 


Fourier space, V • (S a \F tr ue)- We compare it to our es¬ 
timates of the displacement field. As one cannot compute 
the lognormal transform in void regions where 5 = — 1 we 
use the smoothed density field rather than smoothing the 
whole term. This ensures that 6 although converges to never 
reaches -1. 


Scatter plots of the values of V • (S' a SEhrue) versus the 
overdensity and logarithmic transform of the overdensity are 
shown in in Fig. [5] In the higher density samples the lognor¬ 
mal transform is more correlated with the divergence of the 
true displacement field than the 8 field alone. High density 
regions push the 8 values above the values of V • (S + ^true)- 
The lognormal transform suppresses this effect. However in 
void regions where the smoothed overdensity tends towards 
-1 (more so for smaller smoothing scales), the lognormal 
transform varies rapidly to large negative values. 


Fig. m shows a 1 x 50 2 grid point slice through the N- 
body box. The contours show the difference in amplitude 
between V • (S ★ ^true) and both density field estimates 
for three of the samples at both smoothing lengths. The 
improvement from the lognormal transform of 8 compared 
to 8 at estimating the divergence of the true displacement 
field is reduced with sample density. 

In summary, the lognormal transform straightens the 
relationship between density and the divergence of the dis¬ 
placement field for dense samples as previously shown in 


Falck et al. (2012). However, the method suffers from dis¬ 


creteness issues and provides a worse match to the diver¬ 
gence of the true displacement field in low density regions 
compared to the overdensity. These test are only carried out 
on the real space distributions and RSDs are ignored. 

In the second part of this section we test if the lognor¬ 
mal transform of the overdensity provides a better estimate 
of the displacement field, computed using the Zel’dovich ap¬ 
proximation, than the overdensity. We compute the displace- 
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Figure 5. Scatter plots of the amplitude of the divergence of the Lagrangian displacement field at each grid point versus the 5 field and 
the logarithmic transform fields. Different density samples are separated by the horizontal lines, the samples are top to bottom, 20t, It 
and 0p25t. The two plots on the left show the divergence of the Lagrangian displacement field at each grid point versus the the S (far 
left) and versus the the S field and the logarithmic transform (mid left) for smoothing R = 5 h -1 Mpc. The plots on the right show the 
same but for smoothing length R = 10 h ~ 1 Mpc. It is seen that the lognormal transform is better matched to the divergence field in the 
higher density samples. However, in the 0p25t sample when the divergence field is positive, the lognormal transform gets pulled to higher 
amplitude negative values and performs badly. This is more apparent at smaller smoothing scales as expected. Therefore the lognormal 
transform only improves our estimate of V • SI/ in high density regions. 


ment fields from Eq.[2] where both the smoothed overdensity 
and the lognormal transform of the smoothed overdensity 
are used as a proxy for the linear overdensity field. We com¬ 
pare these fields directly to S ★ SFtme- 

Fig. [7] shows contour scatter plots of S * SFtrue * x com¬ 
pared to both \I/true[£ ★ (5(x)] • x an d ^ r true[ln(5' ★ (5(x) + 
1) + C) -x for a subsampled range of smoothing lengths 
and sample densities. We create a new catalogue from the 
N-body simulations that has a realistic galaxy survey aver¬ 
age number density n, where n « 3 x 10 -4 /i 3 Mpc -3 . This 
catalogue has 0.065(128 3 ) particles and is called the 0p065t 
sample. We chose the x direction as the distributions are in 
real space thus isotropic. We quantify our results by com¬ 
puting the standard deviations of the differences between 
these distributions averaged over all cartesian directions, the 
values are shown in Table [4] The displacements computed 
from the smoothed overdensity in the sample with a realis¬ 


tic galaxy survey density, 0p065t, smoothed with a typical 
smoothing length for reconstruction, R = 10 h -1 Mpc, are 
better matched to the true displacements than those com¬ 
puted from the lognormal transform of the smoothed over¬ 
density. 

In Fig. [8] we show the improvement in the estimation of 
^true from a logarithmic transform of the overdensity com¬ 
pared to using the overdensity as a function of sample den¬ 
sity. The values plotted are 1 -ctb/cta, where a a is the stan¬ 
dard deviation of the values of Sl>true — *£(£) and gb is the 
standard deviation of the values of SFtrue — SF(ln(l + 5) + C) 
computed at each grid point. The improvement rapidly 
drops off as the sample becomes more sparse. At the lowest 
density in the 0p065t catalogue, the displacements computed 
from the logarithmic transform of the overdensity smoothed 
with both R — 5 h _1 Mpc and R — 10 h~ x Mpc do not match 
the true smoothed displacements as well as those computed 
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Figure 6. The plots show a 1 X 50 2 grid point slice through the box for the divergence field measured from the 20t density sample 
plus the overdensity field or the lognormal transform of the overdensity for different density samples and smoothing lengths. The best 
match between fields will have the flattest contours and be closest to zero. The 20t sample (top) It sample (middle) and the 0p25t 
sample (bottom) are separated by horizontal lines. The column on the left shows the values smoothed R = 5 h -1 Mpc, and on the right 
smoothed with R = 10 h ~ 1 Mpc. It is seen that for the higher densities, the lognormal transform clearly outperforms the overdensity at 
both smoothing lengths. This can be seen by comparing the hot-spots that emerge in the higher density regions within each sample. As 
sample density is reduced the lognormal transform deviates from the divergence of the displacement field more than the overdensity in 
under dense regions. This can be seen as the blue cold-spots that emerge in the sparse regions of the samples.The lognormal transform 
continues to provide a better match to the divergence of the displacement field in dense regions in the same sample. Reducing the 
smoothing scale increases the differences between the in the two field as expected. 


from the smoothed overdensity. The decline in improvement 
is steeper for smaller smoothing scales as the sample den¬ 
sity is reduced. As realistic galaxy surveys have densities of 
the same magnitude as the 0p065t sample, we conclude that 
it is not beneficial to use the logarithmic transform of the 
overdensity to compute the displacement vectors for recon¬ 
struction. 


7 CONCLUSION 


Reconstruction is an important part of the analysis proce¬ 
dure in galaxy surveys designed to detect the BAO feature. 


Although it is shown to be a robust technique (Burden et al. 


2014), the standard algorithms must be improved for use in 


future galaxy survey data analysis. 

In this paper we have looked at techniques that remove 


the redshift component of the displacement field computed 
from the Zel’dovich approximation applied to the overden¬ 
sity field measured in redshift-space. We have also tested if 
a lognormal transform of the overdensity in the Zel’dovich 
approximation brings the values of the displacement field 
closer to those of the true Lagrangian displacement vectors. 
It should be taken into consideration that tests were carried 
out on a single N-body and mock galaxy catalogue reali¬ 
sation assuming a single cosmology. Furthermore the mock 
galaxy catalogue was created using 2LPT, thus one should 
expect the 1LPT displacement vectors to be recovered to 
higher accuracy than a sample including the full non-linear 
gravitational evolution. We summarise the finding of the pre¬ 
vious sections; 


• The approximation that the redshift component of the 
displacement field is irrotational allows one to compute the 
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Figure 7. Scatter plots of the smoothed Lagrangian displacement field at each grid point versus the displacements computed from the 
smoothed S field and the logarithmic transform of smoothed 5. Different density samples are separated by horizontal lines. The samples 
from top to bottom are 20t, It, 0p25t and 0p065t. The left two plots are smoothed with R = 5 h -1 Mpc and are the displacements 
computed from S and the lognormal transform of S respectively. The plots on the right show the same but for smoothing length 
R = 10 h -1 Mpc. The displacement field computed from the lognormal transform of 5 has a tighter correlation with true displacement 
vectors in the high density samples. However, in the 0p25t sample the displacements calculated from the lognormal transform of the 
overdensity smoothed with R = 5 h -1 Mpc and in the 0p065t sample the displacements calculated from the lognormal transform of the 
overdensity smoothed with R = 5 h~ x Mpc and R= 10 h~ x Mpc do not match the displacement vectors as well as those computed from 
the smoothed overdensity. 
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Figure 8. The improvement in the estimate of the displacement 
field from using a lognormal transform of the overdensity is rep¬ 
resented by (1 — (Jb/&a) X 100 where cr a is the standard de¬ 
viation of the values of ^true — ^(5) and ctb is the standard 
deviation of the values of ^true — ^(ln(l + 5) + C). The im¬ 
provement decreases rapidly with decreasing sample density, h. 
At the lowest density, n ~ 3 X -4 h -1 Mpc -3 , where there are 
0.065(128 3 ) particles in the sample, the displacements computed 
from the logarithmic transform of the overdensity smoothed with 
both R = 10 h -1 Mpc and R = 5h -1 Mpc do not match the 
smoothed true displacements as well as those computed from the 
smoothed overdensity. 


displacement field using FFTs. It is better than no RSD cor¬ 
rection. However, it is shown to overcorrect the RSD slightly. 
Previous work has shown that this does not affect the spheri¬ 


cally averaged power spectrum measurements (Burden et al. 


2014). However it is not clear how reconstruction affects the 


anisotropic statistics. It is possible that the overcorrection 
of RSD due to this approximation may affect quadrupole 
measurements. 

• Empirical computations of the Helmholtz components 
of a plane parallel RSD component to the displacement held 
provide us with an improved approximation. The new ap¬ 
proximation accounts for the irrotational components and 
an estimate of the amplitude of the solenoidal components 
of the RSD term. It is shown to be a better match to the 
real space displacement field in both our N-body simulation 
and mock galaxy catalogue. 

• We propose a fast iterative scheme to recover the real 
space displacement field (provisional on knowing the growth 
rate). Proof of convergence is shown for a plane-parallel RSD 
geometry in Fourier space and we expect the behaviour to 
hold for non-plane parallel geometries. Our iterative method 
does not introduce any distortion in the off-radial compo¬ 
nents in the recovered displacement held as we always start 
from our original overdensity and remove RSD using our 
estimate of \l/ in the radial direction only. The scheme con¬ 
verges in two iterations for the mock sample and N-body 
catalogue when our starting point is the improved approxi¬ 
mation above. The iterative algorithm does not require the 
approximations discussed above, however they speed up con¬ 
vergence. 

• The lognormal transform of the overdenisty helps to 
recover the true displacement held in high density simula¬ 
tions. However is does not perform well when computing the 
displacement held at realistic galaxy survey densities. 


In these tests, the new iterative Fourier algorithm out¬ 
performs the configuration space method of recovering the 


real space displacement field as the accuracy of the configu¬ 
ration space method is limited by the resolution of the grid 
on which the overdensity is computed. However, the config¬ 
uration space method performs very well in terms of the re¬ 
construction process and it has not yet been shown that the 
new Fourier algorithm will offer any significant improvement 
in the effectiveness of the algorithm for current data. The 
effect of the new reconstruction algorithm on the anisotropic 
distribution of the overdenisty has not yet been fully anal¬ 
ysed. It will be interesting to see if the methods described in 
this paper affect the post-reconstruction quadrupole in dif¬ 
ferent ways. We suggest a list of questions to be addressed 
in future work. 

• How well does the iterative FFT reconstruction scheme 
improve the spherically averaged BAO measurement? 

• How do the methods outlined affect the anisotropic 
BAO measurements? 

• Can we use reconstruction to improve RSD measure¬ 
ments? 

As galaxy positions are inferred from their redshift measure¬ 
ments, the background Hubble how and the peculiar mo¬ 
tion of a galaxy due to gravity are coupled along the line of 
sight. This distorts the galaxy clustering map enhancing the 
clustering signal on large scales along the line of sight. The 
anisotropic clustering signal therefore contains information 
about the growth rate of structure. It is possible that the 
reconstruction technique can be manipulated and applied 
to the the RSD held to increase the precision of these mea¬ 
surements. Despite these questions, we have shown that our 
new FFT reconstruction algorithm allows for higher resolu¬ 
tion, more accurate measurements of the displacement held 
from the redshift-space galaxy overdensity than algorithms 
currently in use. This method will be very useful in the anal¬ 
ysis of data expected from the plethora of upcoming galaxy 
surveys. 
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